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ABSTRACT 

A popular formation scenario for giant elliptical galaxies proposes that they might have formed 
from binary mergers of disk galaxies. Difficulties with the scenario that emerged from earlier studies 
included providing the necessary stellar mass and metallicity, maintaining the tight color-magnitude 
relation and avoiding phase space limits. In this paper we revisit the issue and put constraints on the 
binary disc merger scenario based on the stellar populations of disc galaxies. We draw the following 
conclusions: Low redshift coUisionless or gaseous mergers of present day Milky Way like disc galaxies 
do not form present day elliptical galaxies. Binary mergers of the progenitors of present day Milky 
Way like disc galaxies can have evolved into intermediate mass elliptical galaxies (M < M^,) if they 
have merged earlier than « 3-4 Gyrs ago. Assuming that most present day disk galaxies formed in 
a similar way to the Milky Way model presented here, more massive giant ellipticals in general can 
not have formed from binary mergers of the progenitors of present day disc galaxies. A major reason 
for these conclusions is that the mass in metals of typical disk galaxy is approximately a factor of 4-8 
smaller than the mass in metals of a typical early-type galaxy and this ratio grows to larger values 
with increasing redshift. 

Subject headings: galaxy formation: general — galaxy formation: elliptical — methods: analytical 



1. INTRODUCTION 

How did elliptical galaxies form? After more than a 
decade of high resolution observations from space and 
new powerful numerical simulations, additional insight 
has been gained but also new and old questions have been 
raised. Red bulge-dominated galaxies c ontain at least 
half of the stellar mas s of the universe (jFukugita et al.l 
ll998l : lHogg et al.ll2002D and adding bulge s of spiral galax - 
ies the fraction might be as high as 3/4 ( Bell et al]|2003D . 
Therefore a consistent theory for elliptical galaxy forma- 
tion is of fundamental importance but it is still unavail- 
able. Independent of the cosmological context, the lu- 
minous parts of massive galaxies must have formed by 
collecting baryonic material under the influence of grav- 
ity. So, in the most general sense 'mergers' of some type 
are of course necessarily the precursors of present day 
elliptical systems. The questions are when and how they 
have formed. There are a variety of possibilities con- 
cerning how this could have happened. Depending on 
whether the infalling material has already collapsed and 
formed stars the galaxy could either accrete stars or ac- 
crete gas which is then turned into stars thereafter within 
the galaxy. In addition, the distribution of sizes and gas 
fractions of the accreted material might vary with time. 
The choice of these parameters is not completely arbi- 
trary and a detailed theory of elliptical galaxy formation 
must be constrained by comparison with the properties 
of observed elliptical galaxies and of the progenitor com- 
ponents. 

It has been known for decades that giant ellip- 
tical galaxies are a homogeneous fa mily of galaxies 
with old stellar populations (e.g. iFaber fc JacksonI 
IQTetlKormendv fc Diorgovsk|1989l:lThomas et al.l l2005l: 
di Serego Ahghieri et al.ll2006l ). Therefore it has been 
suggested that they all have formed in situ at early 



times (jBrinchmann fc ETIiil200Clf ). Without, at this ear- 
lier time, detailed knowledge of the, currently popular, 
hierarchal cosmological models, the 'monolithic collapse' 
scenario was designed to e xplain the ea r ly for mation of 
ellipticals motivated by the lEggen et al.l (|1962t l proposal 
that the old spheroidal component of the Milky Way 
formed during a short period of quasi-radial collapse of 
gas some 10^*^ Gyrs ago and, in most of their properties, 
the spheroidal bulges of spiral and SO galaxies are indis- 
tinguishable from ellipticals of the same luminosities. 

In that sense, an elliptical/spheroidal galaxy would 
have formed very early, as soon as a finite sufficiently 
over-dense region of gas and dark matter decoupled from 
the expansion of the universe and collapsed. If, dur- 
ing the proto-galactic collapse phase, star formation was 
very efficient, a coeval spheroidal stellar system could 
have formed before the gas could have dissipated its ki- 
netic and potential energy and would have settled into 
the equ atorial plane, thus avoidng the formation o f a disk 
galaxv (iPartridge fc Peebles! Il967t iLarsonl [1961 Il974t 
ISearle et al.lll973l ). Depending on the amount of 'turbu- 
lent viscosity' and the angular momentum of the infalling 
gas, detailed properties like isophotal shapes, rotation, 
age and metallicity gradients and even the formation 
of disk-like substructur es in elliptica l galaxies have been 
computed decades ago (jLarsonl 19751 ) . The 'turbulent vis- 
cosity' was created by large inhomogeneities and random 
motions, eventually caused by independently moving gas 
clouds. Recent high resolution numerical simulations, 
using modern cosmological models, indicate that indeed 
there was a phase of evolution of this type (|Naab et al.l 
|2007( ) during which a small (oc 1 — 2kpc), concentrated 
and massive system was formed rapidly at high redshifts 
(2<z<4). 

At about the same time iToomre fc Toomrd (|1972D : 
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iToomrd ()1974 I1977D investigated the tidal interactions 
during encounters of spiral galaxies and proposed that 
disk galaxies might eventually merge and form ellipti- 
cal galaxies. This proposal became particularly attrac- 
tive with the advent of hierarchical cosmological theories 
(|White fc ReesI (Tots' ) in which mergers play the domi- 
nant role during the formation and/or evolution of every 
dark matter halo and almost every galaxy at some time 
in its life. In this 'merger scenario' massive elliptical 
galaxies could have formed from the gravitational inter- 
action of already existing massive galaxies that were not 
of early type. 

2. THE DISK MERGER HYPOTHESIS 

2.1. CoUisionless disc mergers 

In the early 80's the first fully self consistent simula- 
tions of inerger s of stellar disk galaxies wer e performed 
(iGerhardI [l98il iNegroponte &: White! Il983h^ Using ef- 
fective tree algorithms (iBarnes fc Hutlll98^ iHernquis^ 
[l987t IBarnes fc Hi^ 119891 ^ and" more powerful general 
purpose and special purpose computers a whole indus- 
try of simulations of merging disk galaxies was created 
that was designed to answer the question of whether 
they forr n systems that resemble present day elliptical 
galaxie s jBarncs '1988', '1990; Barnes & Hcrnquist 1992; 
Jl9 



Barnes! [19 92: Hcrn quist 1992: Naab fc Burkcrt 2003 ; 
Jesseit et a l. 2005 :_jGonzalez-Garcia &: Balcellsl 120051 : 
Naab fc Truiillai2006l ). The global properties of the rem- 



nants are in several respects consistent with observations 
of giant elliptical galaxies, e.g. equal mass remnants are 
triaxial, sl owly rotating, an i sotropic, have boxy or discy 
isophotes ((Hevl et all Il99l iNaab fc Burkertll2OO30 . In 
addition, mergers of discs can result in the formation 
of kinematic subsystems like kine matically decoupled 
cores at the center s of e llipticals (jHernquist fc Barnei 
119911 : iJesseit et al.l |2006[ ) as well as observed faint 
structures like shells, loop s and ripples at large radii 
(jHernquist fc Sperge]|ll99"^ . Unequal mass mergers are 
more supported by rotation ([B arnes 1998) and have 
discy isophotes (jNaab fc Burkert) i2003i ) . With respect 
to their intrinsic structure, all coUisionless merger rem- 
nants are dominated by box orbit s at the i r cent ers and 
tube orbits in the outer parts (jBarnesI 119981 ). The 
total fraction of tube orbits increases w ith the mass 
ratio of the mergers (IJesseit et al.l l2005| ). The frac- 
tion of disc and box orbits correlates with the shape 
and kinematics of the systems and the mix of discy 
and boxy isophotal shapes for equal-mass remnants can 
be understood by the pr ojected properties of tube or- 
bits in triaxial potentials ('Jessei t et al.ll200"5t iNaab et al.l 
[20063). [Naab & Burkcrt (2003) and have argued, based 
on statistics of kinematic and photometric properties 
of equal- and unequal-mass mergers, that disc mergers 
(with bulges) can result in intermediate mass elliptical 
galaxies. But the objects so formed are not in agreement 
with t he most massive, bo xy and slowly rotating ellip- 
ticals. iNaab et al.l (|2006bD indicate major binary early 
type mergers could be responsible fo r the slow rotation of 
the rt i ost massive ellipticals (see alsolKhochfar fc BurkertI 
[2005t iGonzalez-Garcfa fc van Albadall200l 

CoUisionless merger remnants in general have phase 
space densities and surface density profiles that resem- 
ble observed ellipticals onl y if bulges are added to the 
progenitor discs (see e.g. ICarlber j Il986f ). The phase 



spac e densities of pure s t ellar disc me rgers are too 
low (jHernquist et all Il993t iNaab fc TruTi llo 2006 ). In 
addit i on, intrinsic and isophotal shapes ((Jesseit et al.l 
120051 : ICox et al.ll2006a|) as well as sur f ace density pro- 
files (jGonzalcz-Garcfa fc BalceUs 20051 iNaab fc Trujillol 
l2006l see however Aceves et al. 2006) of pure disc merg- 
ers are not in agreement with elliptical galaxies. At this 
point, however, there remains the question of how bulges 
have formed in the first place. 

The agreement of the kinematics of coUisionless disc 
mergers and observed elliptical galaxies is only good 
to first order. At higher order there are disagree- 
ments with observed ellipticals. The line-of-sight ve- 
locity distributions (LOSVD) within the effective ra- 
dius of merger remnants in general show small asym- 
metric deviations from Gaussi an shape. They tend 
to have a steep traili i ig wing (iBendo fc Barnes! 12000! : 
Naab fc Burkert! l2001l: iNaab et all l2006al see however 
Gonzalez-Garcia et al. for the influence of massive 



bulges in the progenitor disks), whereas most observed 
rotating e llipticals show a st eep leading wing in their 
LOSVDs ([Bender et al.lll994D . Theoretically, axisym- 
metric, r otating one-component sy stems show such a be- 
haviour (jPehnen fc Gerhardlll994[ ). This would indicate 
that ellipticals are very simple one component systems 
that did not form by mergers, w hich is, however, unlikely 
(see e.g. lEmsellem et al.l [20041 ). An alternative expla- 
nation, based on photometric and kinematical observa- 
tions, is that rotating ellipticals contain embedded large 
scale stellar discs (e.g. Bcndcr'igSft 'Bender ct al."1989!; 
[Rix fc White 1990; Scorza ct al. 1998; Rix et al. 1999). 
A superposition of two distinct components, e.g. a hot 
spheroidal bulge and a rotationally supported cold disc, 
can also result in a steep leading wing of the LOSVD 
([Bender et al.l[l99l INaab fc BurkertI [200lh . 

2.2. Disc mergers with gas 

Evidently, disc galaxies do not only consist of stars 
but also an interstellar medium (ISM) in the form of 
gas. In the local universe for evolved disc galaxies 
typical gas fractio ns are 10 - 30 per cent of the total 
stellar mass (e.g. iMcGaugh fc de Blokl [19971) . Within 
the hierachical paradigm the gas fracti on is an increas- 
ing fu nction with redshift. (see e.g. iKhochfar fc Silk! 
l2006bl ) Even if the overall gas fractions were relatively 
small, due to its dissipative nature gas can change 
the structure of merger remnants sig nificantly (see e.g. 
IKhochfar fc Siia l2006al: iCiotti et al. 12006,) It has been 



shown bv iBarnes fc HernquistI ( 1996 ) that gas accumu- 
lating at the center of merger remnants creates a steep 
cusp in the central potential well resulting in a more ax - 
isymmetric central shape of the remnants (|BarneJ l998') . 
At the same time the fraction of stars on box orbits is sig- 
nificant ly reduced and tubes beco me the dominant orbit 
family ([Barnes fc Hernquist!ll99^ . The most reasonable 
explanation for this behaviour is that systems with steep 
cusps in their p otential can not sustain a l a rge population 
of box orbits (I Gerhard fc Binnev! 1 19851: [Schwarzschildl 
1993; Mcrritt fc Fridman 1996; Vall uri fc MerrittI 119981 : 
Barneail998c iNaab et al.ii2006a) . One interesting con- 



sequence of the change of orbit populations is that the 
asymmetry of the LOSVD of the stars changes in a way 
that is consistent with observations making it a good dy- 
namical tracer for the presence of gas during the merger 
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(|Naab et al.ll2006a[ ). 

Dissipational merging, including star formation, can 
also overcome stellar phase space constraints and there- 
fore the a priori inclusion of a major sheroidal bulge 
component is not required, and it has been shown by 
([Robertson et al.ll2b06a ) that a progenitor gas fraction of 
30 per cent results in remnants parameters in good agree- 
ment with the Fundamental Plane for elliptical galaxies. 
One interesting consequence of gas infall to the center 
and subsequent star formation is a break in the sur- 
face brightness profile with exc ess light at the center 
iMihos fc HernquistlllQQllCox et al. 2006l i an effect that 
is not seen in collisionless mergers ■Naab fc Truiillol[2?M 
For a long time this break in the surface density profiles 
has been considered not to be in agreement with the ma- 
jority of elliptical galaxies. However, recent studies indi- 
cates that a break in the light profile might be a generic 
feature of low and intermediate mass ellitpticals. For ex- 
ample, Sersic functions provide an excellent fit to the pro- 
files of Virgo early-type galaxies outside 10-lOOpc. At the 
centers there is a transition from a luminosity " deficit" to 
a luminosity "excess" with respect to the Ser sic fit mov- 
ing down the early-type lum inosity function (jCote et al.l 
I2006t iFerrarese et al.l l2006lf). It is stil l debated whether 
this transition is smooth ( Cote et al.ll2007[ ) or whether 
it represents a true dichotomy (Kormendy et al. 2007, 
submitted). 

It has been known for some time (|Hernquistl 119891) 
that galaxy mergers can, by dynamical instabilities, 
drive large gas fractions to the center of the rem- 
nants probably causing a starburst and / or feeding 
a super-massive black ho l e (iMihos fc H ernquist 1996; 
Barnes k Hernquisd 119961 : iBournaud e~al.. ,2004 . ,2005, : 
Springel et al.l l2005f ). Processes like this can be ob- 



comparisons. 



3. THE PROBLEM 



served directly in nearby ultra-luminous infrared galaxies 
which are interacti ng disc galaxies with mass ratios in the 
range of 1 :1 to 3:1 (iGenzel et al."1998'. 120011 : iDasvra et all 
[20063) • iHopkins et al.l (|2008. ) have compared surface 
density profiles of observed nearby mergers to simula- 
tions of gas rich disk mergers. They indeed find that 
the observed "excess" central light can be attributed to 
a young stellar component that formed from infalling gas 
during the merger. 

Using a simple model for gas accretion onto a cen- 
tral super-massive black hole it has been argued, based 
on simulations of binary disc mergers, that gas inflow 
regulated by black hole feedback can naturally explain 
the observed present day relation between stellar veloc- 
ity dispersion and black hole mass for ellipt ical galax- 
ies and their evolution with redshift (|Di Matteo et al.l 
120051 [Robertson et al.l l2006b( l. Surprisingly, using the 
same model repeated isolated mergers of more evolved 
systems appear to stay on t he stellar velocity disp ersion 
and black hole mass relation ([Johansson et al.ll2008[ ) . Ex- 
tended models have also been used to understand the 
evolution of Quasar s and stellar spheroids as a whole 
([Hopkins et al.l I2005L [2006) . Howe ver, re cent studies of 
dissipative mergers by (Naab et al] ([2006a) . a nd, includ- 
ing sta r formation and feedback proces ses. bvlCox et al 



2006al lbf) confirm and strengthen the iNaab fc Burkert 



20031 ) conclusion that binary disk mergers only are rea- 



sonable progenitors of intermediate mass giant ellitpicals 
but not of the more massive ellipticals. These studies 
are based on very detailed photometric and kinematic 



Despite the success of the binary merger scenario for 
the formation of elliptical galaxies there remain funda- 
mental problems comparing the stellar populations of 
present day elliptical and spiral galaxies. These prob- 
lems concern even typical ellipticals which would 
from a kinematical point of view be consistent with the 
binary disk merger scenario. Several of the arguments 
pr esented in the fo llowing were outlined in primitive form 
bv lOstrikeii (Il98l . 

One aspect concerns the mass in metals of the two 
kinds of systems. We have estimated the typical mass of 
late-type and early-type galaxies by computing the mass 
above which half the total mass in galaxies of each type 
is contained assuming a lower limit of IO^Mq. Using av- 
erages o f the color and co ncentration selected mass func- 
tions of iBell et all ([20031 ) we get Meariy = 7.5 x 10^°Mq 
and Miate = 2.9 x 10^°Mq. With typical metalhcities for 
the re spective populations o f Zeariy = 0.03 and ziatc — 
0.016 ([Nagamine et al.ll2006( l we get total masses in met- 
als of Mz,early = 2.3 X lO^Mg and Mz,late = 4.6 X IO^Mq. 

E.g. the mass in metals is « 4.8 times higher in typ- 
ical ellipticals than in typical discs. Using an alter- 
native approach using for early ty pe and late type 
galax ies based on the SDSS data of ([Nakamura et al.l 
I2003f ) we obtain M*,eariy = 1-0 x 10"Mq and M*,iate = 

2.5X 10^°MfT, using M/ Learlv = 3.19 and Af/ioarly = 1-18 

([Nagamine eTall 120061 ). With the same typical metal- 
hcities we then get M,n*,z,eariy = 3.1 X IO^Mq and 
M^^^^jate = 4.1 X IQ^Mq. This results in a « 7.7 times 
higher mass in metals of ellipticals. 

Following the model for the evolution of the disc 
of our own ga l axy, the Milky Way, presented in 
iNaab fc OstrikeJ ([2006f ) and extended here we find that 
the stellar component (disc and bulge) has a total mass 
of 5 X lO^^M© containing ft! 8.4 x lO^M© solar masses 
in metals. The mean metallicity is in good agreement 
with typical estimates for the local disk galaxy popula- 
tion as whole, however, the total mass is slightly above 
the typical disc mass estimated before as the Milky Way 
is more massive than a typical disk galaxy. Still we can 
assume that the evolution with redshift is similar for all 
disks. With plausible assumptions for the evolution of 
the star formation rate the mass in metals for one disk 
galaxy was « 50 percent smaller at z = 1 and » 80 per- 
cent smaller at z = 2. Under the simplified assumption 
that early-type galaxies evolve passively since z=2 the 
progenitors of a typical present day disc have about 9-15 
times less metals at z=l and about 20-35 times less at 
z— 2. Simply put the mass in metals of two typical spi- 
rals is significantly less than in typical ellipticals. And 
the problem becomes worse as one considers mergers at 
z=l or even z=2. Clearly this problem is most severe for 
massive ellipticals. 

Another aspect concerns the ages of the stellar pop- 
ulations. There is strong observational evidence that 
evolved, massive, red and metal rich proto-ellipticals 
are already in place at z = 2 — 3 and that present 
day early-type galaxies formed m ost of their stars 

well before a redshi ft z = 1 j Scarle et al. 1973|; 

van Dokkum fc Franxl 119961: iBrinchmann fc Ellis ,200C 



van Dokkum et al.l 119981 : iBender et all 119981 : iTreu et all 
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20051 Ivan der Wei et all 120051: Ivan Dokkum et all 120061: 



Kriek et alJ 120061: iPipino fc Matteuccill2006l ). FoTexam- 
ple. iKriek et al. found, using near-infrared spec- 

troscopy, that about 45% of massive galaxies at z 2 
little/no star formation and evolved underlying stel- 
lar populations. Suprisingly, these galaxies are sig- 
nificantly more compact tha n present day ellipticals 
with equivalen t stellar masses (jvan Dokkum et all 120081 : 
iCimatti et a"n i2008l. 

In general, elliptical galaxies are red and the forma- 
tion of their stars was complete 8-10 Gyrs ago whereas 
disk galaxies in general are blue and have much younger 
stellar populations \ Yith e.g. mean ages of 5 Gyrs 
for the Milky Way (|Robin et al.l 120031) . In combina- 
tion with the fact that disc galaxies are less massive 
and less concentrated than massive elliptical galaxies 
it can be excluded that all elliptical galaxies could 
have formed from binary mergers of present day spi- 
ral galaxies. Of course elliptical galaxies might have 
formed from the progenitors of present day disc galax- 
But observations of the size and mass evolution 



les 



of di s c galaxies (iBarden et al.l 120051 : iTruiillo &: PohlenI 
120051 iTrujillo et al l '20061) as well as a p lethora of 
theoretical models ( Chiap pini et al.l [19971: iMao et al.l 
1998': 'Prantzos & Silk 1998: iBoissier fc P rantzosI l2000l : 
Naab & Ostriker 2006) indicate that spiral galaxies in 
the past were in general smaller and less massive than 
today. As disk merger remnants typically have similar 
projected sizes and not more than two t imes the progen- 
itor mass in the spheroidal component (iNaab fc Trujillol 
[2006) it is even more unlikely that massive elliptical 
galaxies could have formed from binary disc mergers in 
the past. 

Furthermore it has been shown that the stellar popula- 
tions of massive ellipticals have not only formed at high 
redsh i ft but also on s hort ti mescales (e. g. H eavens et al] 
120041 : iThomas et al.l 120051: iPipino fc Matteuccil \200^ . 
This is not compatible with the long formation timescales 
of disc galaxies. The problem is less severe for low 
mass ellipticals which have more extended formation 
timesca les and even show s igns of young stellar popu- 
lations (Thom as et"al] |2005'). 

An additional complication is that elliptical galax- 
ies predomina ntly populate high d ensit y region s like 
galaxy clusters (|Hubble fc Humason|[l93llHubblelll937D . 
whereas disk galaxies populate the field and the over- 
densit y is indepen dent of their luminosity at a given 
color (jHogg et al.l [2003.) . The morphology-density re- 
lationship se ems to be valid over sev e ral magnitudes 
in density ([M elnick fc Sargent I [l977t iD ressW 119801 : 
iPostman fcGel lcr 1984; Whitmo re fc Gilm orc 199^) and 
the most ma ssive elliptic als live in the highest density en- 
vironments ([Hogg et al.li2003) . In summary, the mass of 
spiral galaxies is a very weak function of environment, 
early type galaxies are predominantly found in high den- 
sity regions and the more massive the galaxy the higher 
the overdensity. With the exception of accretion onto 
the central cluster galaxy, mergers are not likely in clus- 
ters, therefore in the simplest form of the merger sce- 
nario a typical elliptical must have formed from 
a merger of two 5 x IO^^Mq disc galaxies outside of the 
cluster, and the hypothesized progenitor spirals are rare 
or nonexistent at any observed epoch. 

In the following we present a model for the evolution 



of a typical spiral galaxy with a central bulge component 
and address several of the above questions using an ide- 
alised scenario of a binary merger of our model galaxy 
and its progenitors. 

4. THE MODEL FOR BULGE AND DISK FORMATION 



In t his section we describe the simple [Naab fc Ostriked 
()2006( l model for the formation of the Milky Way disc 
where we added to our previous treatment the forma- 
tion of a luminous central bulge component. The model 
reproduces nicely observable properties of the present 
day Milky Way without the explicit inclusion of de- 
tailed feedback processes fe.g lEfstathioul l2000l) which 
have t heir main impact at the early phases of disk evo- 
lution Joh ansson fc Efstathiou (2006). We assume that 
in the absence of star formation the gas in a given halo 
would settle in a disc with an exponential surface density 

T,,^d{r,t) = i;o(t)exp[-r/r,nd(t)], (1) 

where the central surface density So(i) and the scale 
length rd{t) can change with time. The galactic disk 
starts to form after the halo has reached its present day 
virial velocity at its formation time im/orm- Thereafter 
it decouples from the general merging. After the decou- 
pling the central surface density is fixed to the present 
day value while the scale length rd{t) is allowed to change 
as a fixed fraction fmr,d of the virial radius rmtiir of the 
halo 



^ md (^ ^ ^ m/orm) — f mr.d 



" (^ m/o 



10H{t) 



(2) 



where 

H[z{t)] = Ho[f7A,o+(l-r!A,o-f^o)(l+-2)'+f^o(l+-2)']'/', 

(3) 

is the Hubble parameter at redshift z. 

At earlier times t < t mform when the galaxy is still 
coupled to the hierarchical growth the gas will have lost 
its angular momentum more effectively due to shocks 
and tidal torques resulting in a smaller collapse frac- 
tion fmr.b- We associate this phase with the formation 
epoch of the bulge. For simplicity we assume that the gas 
also settles in an exponential surface density profile. The 
spheroidal stellar component might form from dynamical 
instabilities or tidal interactions which we do not follow 
explicitly in our model. The scale length of the gas then 
evolves as 

^ mb{i ^ ^ m/orm) — f ni7' mvir (j-^ — f mr.b j^Q^^^^ ^^-^ 

and rotation velocity of the 'bulge' peaks at a radius of 
r m2.2{t) = 2.15rnib(t) at a value of 

<b,2.2W = 0.774^GEo(t)r„b(i). (5) 

Using the similar scalings as in Naa b fc Ostriked ((20M ) 
the central surface density will increase as 

So(t < im/o™) = 4^10i7(i)f ,nc(i), (6) 

with a ~ 3.6 that has been scaled to result in a peak 
rotation velocity the Milky Way bulge of 250fcm/s as- 
suming a present day value of = 210fcm/s at large 
radii. The total surface density of the system has been 
scaled to a present day value of 5OM0pc~^ at r = 8fcpc, 
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Fig. 1. — Upper panel: Time evolution of the total mass, the mass 
in stars and the mass in gas for the model disc. Lower panel: Time 
evolution of the total surface density, the stellar surface density 
and the gas surface density at the solar radius. 
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Fig. 2. — Evolution of the global star formation rate (solid) and 
the global infall rate (dotted). The rapid formation of the bulge at 
3 < 2 < 2 is followed by the quiescent assembly of the disk. 



At this point we have to include a model for star for- 
mation to the model to follow the evolution of the stellar 
and gaseous phase separately. After the gas within a halo 
has settled it starts forming stars. We use a formulation 
based on t he local dynamic al time (rotation period) of 
the system (|Kennicuttlll998D . At every radius the surface 



'mSFR 



with 



{r,t) 



^ ^ ragas (^7 
raorh {Ti ^) 

2-Kr 



(7) 



(8) 



where T,^orb is the orbital period and the star formation 
efficiency is e = 0.1. 

We follow the chemical evolution of the model galaxies 
in independent rings assuming no radial gas flows using 
a modified version of the chemical evolution model of 
lOstriker fc TinslevI (|1975f ). In every ring the change in 
gas surface density S and surface density in stars S 
is given by 



rfS ,„g(r, t) = ~T, ynSFRir, t)dt 
-K^iateir, t)dt, 



^ rains (^i t^dt 

^miFRir, t)dt 
K^ins{r, t)dt 



(9) 



(10) 



where Ti^sfr. is the star formation rate per unit area 
(Eqn. [7]) and Ti^ipn is the rate of gas infall onto the 
galaxy per unit area, as defined in Eqn. [T] i^mins is 
the mass per unit are in gas ejected from massive stars 
instantaneously, -ft'mZate is the mass per unit area in gas 
ejected at later evolutionary phases of low mass stars. 
They are defined as 

(11) 
(12) 
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{r,t)= I S 
Jo 



T,,^SFR{r, t), 

mSFRit',r)Wit^t')dt'. 



R rains — 0.1 is the fraction of gas returned instanta- 
neously to the ISM from newborn massive stars and W{t) 
is a weighting function defined as 

<5. 



W{t) ^ 



1 



To 



To 



t + To 



(13) 



w ith i?» = 0.3, S ^: = 1. 36, and tq = 1 x 1 0^ assuming 
a ISabeteil (fl955l) IMF (jCiotti et al.|[T99l[ ). This ana- 
lytic expression is a good approximation for the frac- 
tion of returned gas for a single burst population to the 
m etal dependent values of th e spectral evolution model 
bv lBruzual fc CharlotI (|2003D that we have used to com- 
pute the photometric propertie s of the m odel disc not 
including the effects of dust (see lNaab fc Ostrikcr 2006). 
We followed the evolution of the metallicity in gas and 
stars both taking into account instantaneous injection 
from massive stars and delayed injection form low mass 
stars. The details of the chemical evolution model are 
given in iNaab fc Ostrikal (|2006l ) . 

4.1. Assembly history 

The upper panel in Fig. [1] shows the evolution of 
the total mass in gas and stars, respectively. The 
galaxy rapidly assembles mass after z = 3, which is the 
bulge formation phase, thereafter the disk is growing. 
At present the total mass of the disc out to 26kpc is 
Mmtot « 6 X lOi°Mm0 with about M^g = 1 x 10^°Mq 
in gas and Mj^tot = 5 x in stars. 

The time evolution of the local surface density is 
shown in the lower plot in Fig. [T] The model has 
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Fig. 3. — Radial metallicity distribution of the gas (solid lines) 
and the mean metallicity of the stars (dotted lines) after 2, 4, 6, 
8, 10, 12 and 13.6 Gyrs (thick lines). The two dotted straight 
lines indicate a slope of —Q.07dexkpc~^ and —0.04:dexkpc~^ . The 
metallicity gradients have been steeper in the past. 



been normalised to a present day total surface density 
of 50M(7)pc~'^ . The gas starts to assemble at the solar 
radius after 3Gyrs, and after GGyrs the local gas surface 
density stays almost constant at its present day value of 
15Mqpc~^ . The first stars at the solar radius form after 
AGyrs followed by a steady increase to the present day 
value of 35Mqpc~^ (with « 3Mqpc~^ invisible in stellar 
remnants) resulting in « 32Mqpc~^ vi sible stars. These 
numb ers are identical to those given in iNaab fc Ostrikerl 
as is the local star formation history at the solar 
radius. Therefore the evolution of the outer disk is not 
affected by the bulge evolution. 

The evolution of the total star formation and infall rate 
of the model galaxy is shown in Fig. [2] The peak in both 
rates at 3 < z < 2 indicates the period of the forma- 
tion of the bulge with infall rates as high as 90MQyr~^ 
and star formation rates of « SOAlQyr^^ . It is interest- 
ing that this analytical treatment of the formation of the 
galactic bulge parallels to a surprising degree the detailed 
hydrodynamical s imulations of an elliptical galaxy by 
iNaab et al.l ()2007l ) from cosmological initial conditions. 
Thereafter the rates drop rapidly to values close to their 
present day numbers of 2 — 3MQyr~^ for the infall rate 
and w 3MQyr~^ for the star formation rate. 

As mentioned before, the evolution of the model galaxy 
at the solar radius is not changed by the inclusion of the 
bulge. Therefore the met allicit y distribution is the same 
and in INaab fc Ostrikerl (|2006l ) and is in good agreement 
with observations. To compute the metallicity evolution 
we have used the solar metallicity value of Zq = 0.0126 
to scale our model to 0.1 dex below the solar metallic- 
ity Z = Zq X 10""-^ = 0.01 at solar radius A.5Gyrs 
ago, resulting in an e ffective yield oi Y ^ 0.0135 (see 
INaab fc Ostri^l200l . 

Fig. |3| shows the radial metallicity distribution of the 
ISM and the stars of the model galaxy. The distribu- 
tion after 2,4,6,8,10, and 12Gyrs is indicated by the 
dotted (gas) and solid (stars) lines. At the solar radius 
th e metallicity grad i ent is similar to the pure disc model 
of lNaab fc Ostrikerl (j2006f ). however, at smaller radii the 
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Fig. 4. — Upper panel: Absolute magnitude of the model galaxy 
versus lookback time. Middle panel: Evolution of the total lumi- 
nosities with time. Bottom panel: Evolution of the color at the solar 
radiu s. The observed colors are B — K = 3.13 ( Binne^_&_^errifielo 
1598) indicated by the star and B — V ~ 0.8 I Boissier fc Prantzoe 
1999i) indicated by the diamond. 
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metallicity of the gas and the stars is significantly higher. 

4.2. Photometric properties 

To compute the p hotometric properti e s of o ur model 
galaxy we use the iBruzual fc CharlotI Hool) models 
for spectral evolution of stellar populations assuming a 
Salpeter IMF. Fig|4| shows the time evolution of the ab- 
solute magnitudes, the total luminosities and the colours 
of the model galaxy at different wavelengths. All val- 
ues are in reasonable agreement with observations. The 
surface brightness profiles at different wavelengths are 
shown in the upper panel of Fig. (5] Two observed values 
in the K- and B-band at th e solar radius are overplotted 
(|Binnev fc Merrifieldiri998l ) . The middle panel shows the 
corresponding luminosity density profiles using the solar 
magnitudes given above. All profiles are exponential in 
the outer parts. The color profiles with two observational 
value at the solar radius are shown in the bottom panel 
of Fig. O We have performed a bulge -disc d ecomposi- 
tion u sing the fitting model described in lNaab fc Truiillol 
(|2006f) and find a bulge to total ratio of B/T=Q.2 in rea- 
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Fig. 5. — Photometric properties versus radius. Observed 
values are HmK = 20.62 (star) and fimB = 23.75 (diamond) 
(iBinnev & Merrifield 1998) The observed colors are B — K = 3.13 
ijBinncv & M errifiel d 1998) indic ated by the star in the lowest 
panel and 0.6 < B - V < 0.8 (see|B oissier fc Prantzo3[T999h indi- 
cated by the error bar. 



sonable agreement with estimates for the Milky Way. 

To calculate the scale lengths in the different bands we 
have fitted an exponential in the range of 1.5 < r^^^d < 
3.0 scale lengths of the total surface density distribution. 
The scale lengths increase to shorter wavelengths (Fig. 
[H]). This effect is weaker at earlier times. At present the 
B-band scale length r ^d,B = 3-7kpc is a factor of « 1.2 
larger than the K-band scale length r^d,K = 3.0A:pc . 
This trend is observed and is in good agreement with 
other Milky Way type spiral galaxies. Ide Jongj ()1996f ) 
investigated a sample of 86 nearly face-on spiral galax- 
ies and concluded that spiral galaxies of the same type 
as the Milky Way have disc scale lengths that are a 
factor of 1.25 ± 0.25 larger in the B-band than in the 
K-band. This effect reflects the over all inside out for- 
mation process for disc galaxies and has b een found in 
similar studies (jBoissier fc: Prantzoi (|1999[ ) have found 
r md.s/r nid,K = 1-5 which is also in good agreement). 
We have separately examined the age distribution of the 
stars in the model Milky Way galaxy and find them in 
good agreement with observations, giving additional cre- 
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Fig. 6. — Exponential scale lengths of the surface brightness 
distribution measured at the different wavelengths versus time. At 
shorter wavelengths the disc has a larger scale length. 

dence to our believes that the backward evolution of the 
model galaxy is plausible. To summarize the evolution 
of our model Milky Way, it was smaller and less massive 
in the past but since z 2 its luminosity and metallicity 
have been remarkably constant. 

5. THE IDEALISED DISC MERGER 

Numerical studies in the past have shown that sizes, 
photometric properties and dynamics of disk merger rem- 
nants are in reasonable agreement with the properties 
of observed intermediate mass giant elliptical galaxies. 
Using the model for galaxy evolution described in the 
previous section we now have a self consistent model 
giving knowledge about the evolution of a typical disc 
galaxy, e.g. ages and metallicities of the stars and gas 
content at any redshift. Using this information we can 
investigate the remnant properties in an idealised disk 
galaxy merger scenario. We assume that two identical 
disk galaxies (with or without bulges) merge at different 
times in the past and investigate the integrated proper- 
ties of their combined stellar populations today. During 
the merger we allow all gas within two disc scale lengths 
to be transformed into stars in a single burst. The dura- 
tion of the burst was taken from typical burst timescales 
of equal mass binary disk merger simulations (see e.g. 
[Johansson et al.l [20081 ). Somewhat smaller or larger du- 
rations do not influence the results. After the galaxies 
have merged we assume that further star formation is 
suppressed. This might either be due to the heating of 
gas by supernovae, AGN feedback, shock heating of in- 
falling gas, the infall of the remnant into a cluster or a 
combination of all effects. As we know the properties of 
the stars and the gas content of the progenitors disks at 
any redshift, we can predict the present day properties 
of the stellar population of the remnant that has formed 
by a merger at any time in the past. 

In Figure [7] we show the present day location of the 
merger remnants of two Milky Way like galaxies in the 
U — V color-magnitude diagram. For every merger model 
the rightmost symbols represents the properties of the 
present day merger. Following the symbols for every 
model from the right to the left we have placed the 
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Fig. 7. — Upper panel: U — V color versus absolute magnitude 
Aiv of the stellar populations of two pure disk galaxies added to- 
gether in an idealized merger at different times in the past. The 
triangles show the properties for the stellar remnant alone. The as- 
terixes show the remnant properties assuming that all gas within 
two scale lengths at the given time was transformed into stars in a 
single burst with a duration of 150 Myrs. The burst and no-burst 
models are connected by a dotted line. The present day burst 
model is off the scale. The rightmost pair indicates the properties 
of a present day merger remnant. Every pair further to the left 
shifts the merger event one Gyr further back into the past. Lower 
panel: Same as above but for for progenitor disks with bulges. All 
remnants are consistent with the CM-relation if the merger would 
have taken place 3 Gyrs ago or earlier. 5 Gyrs or older mergers of 
disk galaxies with bulges lead to the most consistent results. Merg- 
ers of present day disks do not fall on the CM-relations independent 
of their properties. 



merger further back in the past. The observed UV color- 
maRnitude diagram o f eUiptical galaxies is overplotted 
(|McIntosh et al.ll2005D . In the two upper panels we show 
the properties of a me rger of pure disks with and without 
a burst assuming the iNaab fc Ostrikerl ((2006) model. It 
is evident that a present day merger of disk galaxies re- 
sults in a remnant that is far too blue to be an elliptical 
galaxy. The discrepancy is even stronger if we consider 
the young population of stars created in a burst dur- 
ing the merger. Placing the merger further back in time 
results in a present day remnant that becomes less lumi- 
nous and redder. Taking the spread in the CM-relation 
into account the remnant properties would be consis- 
tent with both CM-rclations if the galaxies would have 
merged 3 — 8 Gyrs ago. The remnant would fall right on 
the CM-relations if the merger would have taken place 
5 Gyrs ago and would give rise to an elliptical galaxy 
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Fig. 8. — Mean age of the stellar population of disk mergers 
with and without bulges versus the lookback time of the merger. 
The horizontal line indicates a lower limit of 8 Gyrs for the age 
of giant elliptical galaxies. All younger ages (shaded area) are not 
compatible with the stellar populations of elliptical galaxies. 



with an absolute V-band magnitude of A/y — —19.8 at 
present. Surprisingly, if we go even further back in time 
the remnant of the model without a bulge would lie above 
U — V CM-relation leaving two possible interpretations: 
either the remnant is too red or it is simply not luminous 
enough. 

If we add a bulge to the progenitor disk using the model 
presented in the previous section the remnants become 
redder and more luminous at all times (lower panel in 
Fig. [T]). As in the disk-only case we have to place the 
merger event at least 3 Gyrs in the past to be consistent 
with the present day UV CM-relation. The 5 Gyr old 
merger falls right on the correlation as in the pure disk 
case and produces an elliptical having a current mag- 
nitude of My = —20.6. In contrast to the disk-only 
models, all earlier mergers fall on the color-magnitude 
relations due to the presence of the massive old bulge 
component. However, at higher redshifts 2; > 1 it is a 
merger of two bulge dominated galaxies and not a disc 
merger, i.e. to a significant extent we are considering the 
merger of ellipticals to make new elliptical. 

In Fig. [8] we show the mean age of the stellar popu- 
lation of the merger remnant at different lookback times 
for the models with and without a bulge component. For 
a 5 Gyr old merger (which corresponds to a merger red- 
shift of z « 0.5) the disk-only remnant has a mean age 
of 7.9 Gyrs. The remnant with bulges in the progenitor 
disks is significantly older at a mean age of 9.2 Gyrs. 

Our model Milky Way galaxy with bulge is about 
two times more massive than a typical present day disc 
galaxy. If wc, however, assume that the backward evo- 
lution with time is typical for disk galaxies of simi- 
lar mass, and there is observational evidence for that 
(|Barden et all 120051 ). we can scale our model down and 
follow the evolution of the mass in metals for a typi- 
cal disk merger backwards in time. In Fig. [9] we have 
scaled the model to a present day total stellar mass of 
2.9 X IO^^Mq (see Section [3|) and compare the mass in 
metals of disc mergers with and without bulge to present 
day ellipticals. At present, disc merger remnants have 
« 2.4 — 3.1 times less mass in metals if they contain 
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Fig. 9. — Evolution of the mass in metals for merger of two 
Mt disk galaxies based on the model presented in the paper with 
(solid) and without (dashed) a bulge component. The total mass 
of the progenitor disks has been scaled to 2.9 X IQ^^Mq. The mass 
in metals for Af, ellipticals is indicated by the shaded area (see 
Sec. [Sjl- At present ellipticals have at least a factor of two more 
mass in metals. 

bulge components. Without bulges this value becomes 
as large as 6-8 at the present day and even larger at 
higher redshifts. The bulge, although less massive than 
the disk itself is significantly more metal rich, in line 
with our previous arguments. Therefore the total mass 
in metals is about a factor of two higher for the model 
with bulge. Assuming a burst during the merger does not 
change the overall values much as the total gas fraction 
of the progenitor galaxies stays at a constant fraction of 
« 20% since z—2. The metal mass of ellipticals could 
only be reached if all available gas would be transformed 
into stars with eight times the solar metallicity during 
the merger. 

6. CONCLUSION AND DISCUSSION 

In this paper we have addressed the question of 
whether elliptical galaxies can have formed from binary 
mergers of disc galaxies. At present typical disc galax- 
ies have lower stellar masses, lower masses in metals, 
younger stellar populations and more extended star for- 
mation histories than elliptical galaxies. In addition they 
populate different environments. The first simple con- 
clusion is that mergers of present day spiral galaxies, 
like our Milky Way, can not form typical present day 
ellipticals but some subset of future ellipticals having 
L < may form in this fashion right now. In fact 
nearb y ULIRGS are disc m ergers with mass ratios of 1:1 
to 3:1 (jPasvra et al.ll2006al ) and they show kinematic and 
isophotal properties similar to present day intermediate 
mass ellipticals ((Ge nzcl ct al. 2001; Rothbcrg & Joscg. 
120041: iDasvra et aH I2006al lbl: IRothberg fc JosephI |200" ° 
After their stellar populations have aged for r:! 3 — AGyrs 
they will add to the population of intermediate mass el- 



lipticals. Expanding the investigation to cosmic history 
we have used a theoretical model for mergers of a typical 
Milky Way spiral at different epochs to put more general 
constraints on binary merger scenario. 

Assuming that most present day spirals formed in way 
similar to our model presented here we can place further 
constraints on the formation mechanisms of early-type 
galaxies. 

We find that mergers of progenitors of present day spi- 
ral disc galaxies can have stellar populations compara- 
ble to present day low and intermediate mass elliptical 
galaxies if they have merged more than 3-4 Gyrs ago 
(z > 0.3) and have their star formation turned off after 
the merger event. Any further star formation after the 
merger will shift the time for the merger event further 
into the past. The progenitors of mergers before z=0.3, 
however, have smaller sizes (r < 2.5fcpc), smaller stellar 
masses (M < 4x IO^^Mq) and are more bulge dominated 
than at present. As we know from a large number of nu- 
merical studies (see Section these will have sizes and 
kinematical properties similar to intermediate mass (M*) 
ellipticals. And early mergers of spiral galaxies will lead 
to ellipticals having metal masses far below those seen in 
present day ellipticals. 

Specifically, if the merger takes place further back in 
the past, say at z=l, the typical sizes of discs will be 
(r < 2kpc) with stellar masses M < 3 x 10^°Mq. How- 
ever, the typical mass in metals will be a factor of 4-8 
smaller than typical present day ellipticals, and, even 
worse, at z = 2 the mass in metals will be a factor of 
10-18 smaller and the merger progenitors will be small 
and bulge dominated. This epoch, however, is supposed 
to be the main formation epoch of ellipticals, especially 
the most massive ones. 

Again, assuming that most present day disks formed in 
a similar way, we can conclude that typical massive gi- 
ant ellipticals (more massive than Ms,) can neither be 
made from binary mergers of present day spirals nor 
from their progenitors although some low mass subset 
(less massive than M*) of the total observed z=0 ellip- 
tical galaxy population is certainly produced by such 
mergers. Alternatively, massive ellipticals coul d have 
formed by either multiple mergers (|Li et al.ll200"l or bi- 
nary mergers of massive spirals at early times whose de- 
scendents no longer exist. These massive high redshift 
spirals would have had different formation histories than 
the model galaxy presented here. Any scenario address- 
ing the full mass range of ellipticals that is based on 
binary mergers, e.g. for the origi n of the black-hole-mass 
velocity-dispersion rel ation (e.g. iDi Matteo et al.l [20051 . 
[Johansson et al.ll2008[ ). will have to account for this con- 
clusion. However, there is clear evidence that dissipation 
played an important role during the forma tion of the 
spher o ids especially at low e r masses (see e.g. iKormendyl 



19891: iBender et a.1.1 119921: iKormendy fc Bender T99 



n n n H- 

Dekel fc Coxl 120061: Naab et all l2006air iKh ochfar & Sif 

2006ai iCiotti et al.' 20061: iNaab et al.l [2007)1 . However , 
even massive ellipticals with rotation supported, old and 
metal rich kine matically decoupled cores show clear signs 
of dissipation ([Bend er fc Surma' 199^; 'S urm a fc Bended 
ll995tlDavies et al.l l2001: McDcrmid et al. 200^^ 

Major binary mergers of early- type galaxies might play 
an i mportant role for the final assembly of massive ellipti- 
cals (jkhochfar fc Burkertj[2003t [Robertson et al.ll2006al : 



Naab et al. 



10 

Boylan- Kolchin et al.ll2006HBell et al.ll2006al:lNaab et all 
2006b; Be ll et al]l2006bl ) but they cannot be connected 
to the formation of the spheroids in the first place. 
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